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I. INTRODUCTION 



In this paper we demonstrate the feasibility of QCD sum rule calculations |I[ for CP-odd 
electromagnetic observables induced by the QCD vacuum angle 9. This parameter labels 
different super-selection sectors for the QCD vacuum, and enters in front of the additional 
term in the QCD Lagrangian 

g 2 

(-> = 9 rG" ,G .... (1) 

3271-2 ^ v w 

which violates P and CP symmetries. As it is a total derivative, G a G a can induce physical 
observables only through non-perturbative effects. 

Experimental tests of CP symmetry suggest that the ^-parameter is small, and among 
different CP-violating observables, the electric dipole moment (EDM) of the neutron is one 
of the most sensitive to the value of 9 |J. The calculation of the neutron EDM induced by 
the theta term is a long standing problem. According to Ref. |§, an estimate of d n (6) can 
be obtained within chiral perturbation theory. The result, 

, a m u m d ,0.9 

dn = - e6 mm u + md ) { 4^ + C) (2) 

is seemingly justified near the chiral limit where the logarithmic term becomes large and may 
dominate over other possible contributions parametrized in this formula by the constant c. 
This constant is not calculable within this formalism and in principle can be more important 
numerically than the logarithmic piece away from the chiral limit. In fact it is also worth 
noting that in the limit m u , — > 0, the logarithm is still finite, and stabilized, for example, 
by the electromagnetic mass difference between the proton and neutron. In any case, an 
inability to determine the size of corrections to the logarithm means that one is unable to 
estimate the uncertainty of this prediction. 

If the logarithm is cut off at A ~ m p and non- logarithmic terms are ignored, one can 
derive the following bound on the value of 9 using current experimental results on the EDM 
of the neutron ||] : 

#<3-l(T 10 . (3) 

Confronted with a naive expectation 9 ~ 1, the experimental evidence for a small if not zero 
value for 9 constitutes a serious fine tuning problem, usually referred to as the strong CP 
problem. The most popular solution for the strong CP problem is to allow the dynamical 
adjustment of 9 to zero through the axion mechanism |J. 

There are two main motivations for improving the calculation of the EDM of the neutron 
induced by the theta term. The first refers to theories where the axion mechanism is absent 
and the ^-parameter is zero at tree level as a result of exact P or CP symmetries |||7|]. At 
a certain mass scale these symmetries are spontaneously broken and a nonzero 9 is induced 
through radiative corrections. At low energies, a radiatively induced theta term is the main 
source for the EDM of the neutron as other, higher dimensional, operators are negligibly 
small. As 9 itself can be reliably calculated when the model is specified, the main uncertainty 
in predicting the EDM comes from the calculation of d n (9). 
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The second, and perhaps the dominant, incentive to refine the calculation of d n (6) is 
due to efforts to limit CP-violating phases in supersymmetric theories in general, and in the 
Minimal Supersymmetric Standard Model (MSSM) in particular. Substantial CP-violating 
SUSY phases contribute significantly to 9 and therefore these models apparently require the 
existence of the axion mechanism. However, this does not mean that the ^-parameter is 
identically zero. While removing 9 ~ 1, the axion vacuum will adjust itself to the minimum 
dictated by the presence of higher dimensional CP-violating operators which generate terms 
in the axionic potential linear in 9. This induced ^-parameter is then given by: 

9 ind uced = -^t, where K x = i { / dxe^ (0\T(^GG(x), O CP (0)\0)\ , 

\K\ U 671 J q=0 

where Ocp{0) can be any CP- violating operator with dim>4 composed from quark and 
gluon fields, while 

K = i { / dxe^(0\T(^GG(x), ^GG(0))|0)} (4) 

is the topological susceptibility correlator. In the case of the MSSM, the most important 
operators of this kind are colour electric dipole moments of light quarks qig^G^a^^q, 
and three-gluon CP-violating operators. The topological susceptibility correlator K was 
calculated in f|||J and the value of 9, generated by color EDMs can be found in a similar 



way flOfl . Numerically, the contribution to the neutron EDM, induced by 9 e jf is of the 
same order as direct contributions mediated by these operators and by the EDMs of quarks. 
Therefore, the complete calculation of d n as a function of the SUSY CP-violating phases must 
include a d n (9) contribution and a computation of this value, beyond the naive logarithmic 
estimate (Q), is needed. 

Within the currently available techniques for the study of hadronic physics, it seems 
that the only chance to improve analytically on the estimate @ is by use of the QCD sum 
rule method Given its success in predicting various hadronic properties, including the 
electromagnetic form factors of baryons > ft appears highly suitable for the calculation 

of observables depending on 9. In the sum rule approach, physical properties of the hadronic 
resonances can be expressed through a combination of perturbative and nonperturbative 
contributions, the latter parametrized in terms of vacuum quark-gluon condensates. In the 
case of CP-odd observables induced by 9, the purely perturbative piece is absent and the 
result must be reducible to a set of the vacuum condensates taken in the electromagnetic and 
"topological" background. The expansion to first order in 9 will result in the appearance of 
correlators which have a structure similar to K and K\ in Eq. ([|). These correlators can 
then be calculated via the use of current algebra, in a similar manner to those considered 
in HH]. In this approach the ^-dependence will arise naturally, with the correct quark 
mass dependence, and the relation to the U(l) problem will be explicit. This relation is 
manifest in the vanishing of any GG-induced observable in the limit when the mass of the 
U(l) "Goldstone boson" is set equal to the mass of pion. 

Obviously, the calculation of the EDM of the neutron induced by the theta term will 
be a substantial task, although it appears that the main problem may be technical rather 
than conceptual - the calculation of d n (9) needs to be at linear order in the quark mass as 
compared to the calculation of the anomalous magnetic moment which may be performed in 
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the limit m U) d = 0. There are, however, additional subtleties relating to the correct choice 
for the nucleon current in the presence of non-zero 9. 

Keeping in mind the importance of a sum rule calculation for the EDM of the neutron, 
we would like to test the applicability of this method by calculating the EDM(#) for a sim- 
pler system. The perfect candidate for this would be p meson which couples to the isovector 
vector current and whose properties have been predicted within QCD sum rules with im- 
pressive accuracy [TJ]. Thus in this paper we propose to study the feasibility of sum rule 
calculations for ^-dependent electromagnetic observables in this mesonic system. We begin, 
in Section 2, with a tree-level analysis of the vector current correlator in a background with 
a nonzero electromagnetic field, and a theta term. We obtain the Wilson OPE coefficients 
for all ^-dependent contributions up to neglected operators whose momentum dependence 
is 0(l/q G ). This result is also briefly contrasted with the analogous expression for the tensor 
structure leading to the magnetic dipole moment. In Section 3, we turn to the phenomeno- 
logical side of the Borel sum rule, and in Section 4 study the various contributions to the 
sum rule in some detail, obtaining a stable relation at the level of 0(l/g 4 ) terms, and a rea- 
sonably precise extraction of d p . Section 5 contains further discussion, including comments 
on a consistent procedure for the definition of current operators away from the chiral limit 
in a background with nonzero 9. 



II. CALCULATION OF THE WILSON OPE COEFFICIENTS 

Since it is a spin 1 particle, the p-meson can possess on-shell two CP-odd electromagnetic 
form factors, the electric dipole and magnetic quadrupole moments. We shall concentrate 
here on the EDM of p +( -~^ as the CP- violating form factors of p°, induced by the theta term 
should vanish. This is a general consequence of C-symmetry, which is respected by the theta 
term. 

Before commencing the calculation, it is useful to have a rough estimate of the EDM of 
p induced by theta. It is clear that the correct answer should have the "built-in" feature of 
vanishing when m u or rrid are sent to zero. Thus we should expect a result of the form, 

e m u m d 

dp ~ — — -, (5) 

m p A(m u + m d ) 

where A is some scale at which the reduced quark mass is effectively normalized, presumably 
between f n and m p . We note in passing that one could also use the approach of [[| to obtain 
the contribution to d p due to the chiral logarithm. 

In order to calculate the p + EDM within the sum rule approach, we need to consider the 
correlator of currents with p + quantum numbers, in a background with nonzero 9 and an 
electromagnetic field F pu , 

IWQ 2 ) = i J rf 4 xe^(0|T{j+(x)j;(0)}|0),, F , (6) 

where we denote Q 2 = —q 2 , with q the current momentum. 

To simplify the presentation, we shall consider the example of p + in the m u = limit, 
for which the current reduces to j+ = uj^d. Since we always work to linear order in the 
quark mass, it is straightforward to resurrect the full dependence when required below, and 
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we shall always write the full mass dependence explicitly, with the implicit understanding 
that we set m u = m<j. With this current structure, the correlator reduces to the form 

n+,(g 2 ) =i J dW^(0|n(x) 7 ^ d (x, 0)7^(0) + d(0)^ u S u (0,x) llx d(x)\0)g, F 

= + rfu,, ( 7 ) 

where we have contracted two of the quark lines leading to the presence of the d and u quark 
propagators, S u (x, 0) and S d (x,0), respectively 

We concentrate now on the contribution 7r" , and note that at tree level the linear 
dependence on the background field can arise either through a vacuum condensate, 
or from a vertex with the propagator. These contributions are depicted in Fig. 1, and 
correspond to an expansion of the quark propagator to linear order in the background field. 
If we assume a constant field dpF^ = 0, the gauge potential may be written in covariant 
form Afj,(x) = — ^Fp U (0)x u , while similarly if we work in a fixed point gauge fll3| , the gluon 



gauge potential may also be represented as A^i^x) = — |G^(0)t a x^. The expansion of the 
massless propagator, conveniently written in momentum space, then takes the form W4 



S(q) = J d'xe^Six, 0) = - + ^- e P am5 + ^- g G a p W + ■ • • , (8) 

where G^ v = G^t" 1 , and we have introduced the dual field strengths F a p(— ^e a p pcT F pCT ), and 
Gap. Since we are concerned only with the leading linear dependence on the quark mass, the 
mass structure of the propagator is very simple, and while not shown explicitly here, this 
structure is easily resurrected when required. The particular contributions we shall need 
will be given below. 

While the expansion of the propagator (H) apparently exhausts all possibilities for ob- 
taining a linear dependence on the background field, it is important to also consider an 
expansion of the quark wavef unctions. The first order correction in the covariant Taylor 
expansion will be sufficient here, and is given by 

u(x) = tt(0) + x a D a u(0) H , (9) 

where D a = d a — ie u A a (x) is the covariant derivative in the background fieldP]. 

Thus, if we substitute the expansions for the quark wavefunction and propagator into 
7r" , we find a sum of six terms, 

C = * K + 4 + 4 + + « u 82 + , (io) 

which may be conveniently represented in momentum space as follows: The first three terms, 



1 In principle, next-to-leading order contributions linearly proportional to the field strength may 
also arise from the second order term in the Taylor expansion. However, these contributions have 
a very small coefficient [11|, and will be ignored, as they constitute a negligible correction to the 
terms we shall discuss below. 
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FIGURES 




< = (0|«(0) 7 ^7X0)|(W (11) 

tt 2 u = g(0\u(0)j,^-(Ga) lv u(0)\0} e , F (12) 

irt = eMmi,^ 2 (Fa) lu u{m)e (13) 

represent the contributions at leading order in the quark wavefunction, while ~Kq X 2 3 corre- 
spond to the first order corrections, 

= -i<0|n(0) D a 7 , - 7,7^(0)10)^ (14) 

tt£ 2 = -fc(0|B(0) D a 7/1 (tJ^ " G^(0) 77 767,«(0)|0>^ (15) 

7& = -ie d (0\u(0) D a lfl ( j% - ^,(0)7^(0) |0) F . (16) 



Note that in evaluation of 7r 2 and 7Tg 2 we may turn off the background electromagnetic field. 

We shall now consider each of these contributions in turn, although as a first step its 
helpful to re-express the quark wavefunction corrections 7Tg x 2 3 in terms of the leading order 
condensates via use of the equations of motion. It proves convenient to first analyze the 
derivative term 71"^. Writing this result in spinor notation, so as to factorise the 7-matrix 

structure, the matrix element we need to consider has the formQ (0|w a D a Ub\0). Since we are 
only concerned with matrix elements which give a nonzero contribution when evaluated in 
the ^-vacuum, its helpful to choose an appropriate basis of condensates in which to expand. 
For this example, there are two natural vector and axial-vector structures to consider. We 
write 

(0\lI a D a U b \0) = Ci(7 /3 )6a(0|Me a/ 3 M A-D At 7A75w|0) 

+C 2 (lf3l5)ba(0\uD [a >y l3] >y 5 u\0), (17) 

where C\ and C 2 are two constants to be determined and [a, (3] is used to denote anti- 
symmetrisation of indices (symmetrisation of indices will later be denoted by {a, /?}). 



2 In principle one could also explicitly include colour indices. However, the trace over these indices 
will always be trivial in the examples to be considered here, so this dependence will be suppressed. 
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Making use of the following identities, 

{^,^} = -2e a ^V75^ (18) 
W>, v*f3] 75 = 2iD [a 7 /3] 75, (19) 

integrating by parts, and using the Dirac equation ]J)u = —im u u, we can reduce the natural 
Lorentz decomposition (|T7j) to a more recognizable form, 



(0\u a D a u b \0) = im u Ci(7 /3 ); )a (0|M(j a/ 3M|0) 

-m u C 2 ('jp'J5)ba(0\ua a( 3'y 5 u\0). (20) 

Two equations for the constants C 2 and C 2 may be obtained by, in one case, contracting 
with (a^ajab, and in another, via multiplication by 75 and then anti-symmetrising in a, 5. 
One obtains C\ — and C 2 = —1/8, and thus we are left with only one structure in the 
decomposition. 

On Fourier transforming to momentum space, and performing the straightforward 7- 
matrix algebra we find 



7T 



-= z^(0|uo-^w|0) eii? - m d ^Je^ 7/3 (0|M(T Q/ 375M|0) eiF . (21) 



We can now compare the first term here with 7r" given in (|TT|) . Since the ^-dependent CP- 
odd contribution corresponds to considering only the term antisymmetric in and v, the 
relation 7^7^ = — icr^+sym. implies that 7r" precisely cancels the first term above. Thus 
we have, 



7Ti + 7T. 



m d p^re flull3 (0\ua a ^ 5 u\0) eiF . (22) 



Turning next to ttq 2 (|T^), a little 7-matrix algebra shows that the matrix element may 
be rewritten in the form 

= (j^y 2 ~ e^ x F pa (0\u D a lx u\0) e . (23) 

It can be shown that the condensate in this expression is in fact proportional only to 
5 , aA^u(0|H'u|0)6i and therefore does not contain any CP-violating piece. 

The final wavefunction correction to consider is ttq 3 (16) which may be handled in a 



similar manner to ir^, via extracting the appropriate projections onto vacuum condensates, 
or alternatively by direct calculation. We shall follow the former approach here, and write 
down 

(0\u p (G a ^)D a u q \0) = C x {t a ){ la ) pq {U\u{G M ) #m|0). (24) 

Note that another apparently valid Lorentz structure of the form (0|u(G^ 7 )7 5 flu\0) vanishes 
on the equations of motion. Contracting with (7 a )<j P (£ a ).,-j, and recalling that t a t a = 4/3, one 
finds that C\ = 3/64. The resulting expression for 7ig 3 is, 



n 03 ' A I ^2 ^ 2V! 



4= e^ a7 (0|uG^|0). (25) 



The next problem to address is that of extracting the leading ^-dependence of these 
matrix elements, and we follow standard practice (see e.g. ||) in making use of the anomalous 
Ward identity. To illustrate the procedure, consider the condensate m u (0|wrw|0), with a 
generic Lorentz structure denoted by T. In the vacuum, to leading order, we have 

m u (0\uTu\0) e = m u J d 4 y(0\T{(uTu(0), i9^-G%G a ^ (y)}\0) . (26) 

We now make use of the anomalous Ward identity for axial currents || restricted to 2 
flavours. A useful calculational simplification follows if we take as the anomaly relation a 
linear combination of the singlet equations for the u and d quarks. In particular, we use 



where 



= 2m*(ui l5 u + dt l5 d) + ^-G%G a »\ (27) 



= —ui^u + —dj^ 5 d (28) 
m u m d 



is the anomalous current, and we have introduced the reduced mass, 



m* = (29) 

m u + m d 



Substituting the anomaly relation ( P7|) for GG into the correlator (|26[ ), we recall that 
the only contribution from d^j^ is a contact term oc S(y ) due to the presence of the T- 
product. Through the use of the equal time commutator, we find that this leads to a local 
contribution (independent of y), and consequently we have 

Tfl 

m(0\uTu\0)e = iOm u —{Q\uT^u\$) 

+16 J d 4 3/(0|r{m,(«7 5 tt(j/) +d l5 d(y)),m u ua„ u u(0)}\0). (30) 

The nonlocal contribution to this correlator, the second term above, is 0(m 2 ) in light quark 
masses. Nonetheless, this term would cancel the local contribution were there an interme- 
diate state with mass squared of 0(m) - for example the Goldstone boson in the singlet 
channel. The crucial point, as stressed in |§, is that due to the U(l) problem the lightest in- 
termediate state r] has m n 3> m„- and thus the second term can be neglected at leading order 
in m. Thus for each of the contributions above, the leading dependence on 9 is determined 
via the following relation, 

m(O\uTu\O) = i9m^(0\uT^ 5 u\0) . (31) 

We could of course have obtained this result in a simple manner by using the anomaly to 
rotate away the GG term in the action. This induces a complex quark mass m — > m+iOm^s, 



and leads directly to the leading ^-dependence ([!!]) above. However, despite being somewhat 
more involved, the procedure we have followed is advantageous in that it makes quite explicit 
the role of the anomaly and, in particular, the conditions under which the higher order non- 
local terms may be neglected. We shall return to the issue of chiral rotations at the level of 
the action in Section 5. 



S 



The final effect to consider is that of the background field, F f 



For the term ir^ 



OD, the 



leading F-dependence is already explicit, and may be extracted via introduction of spinor 
indices. For the other terms, we follow Ioffe and Smilga 
susceptibilities", Xi K > an d £, defined as follows [|Tl"11 : 



11 and introduce "condensate 



(0|g^g|0) F = e qX F^(0\qq\0) 
9(0\q(G; u t a )q\0) F = e q KF, u (0\qq\0) 
ge^ a (0\q l5 (G a X(T t a )q\0) F = ie q £F, u (0\qq\0), 



(32) 



where q = u or d. 

Using the relations 



llj32|), and performing the Fourier transformation to momentum 



space, we can now gather all the results from vr"-7rg 3 
analogous results for k a in (|7|), to obtain 



TO!) and combine them with the 



11+ =m*0(e tt -e d )<O|gg|O) 



hi/ 



-x 




- x 



Q 



2\2 



(33) 



This expression exhibits the two tensor structures one would have expected to appear on 
general grounds. However, only the first term contributes to the EDM, as one may check by 
choosing a rest frame for the current momentum, since the second tensor structure vanishes 
on shell. 

Therefore, if we retain only the contribution which survives on-shell, we have as our final 
expression for the theoretical side of sum rule, 



n 



fJ,V 



mj(e u - e d )(0\qq\0) 



q 



-x 



l 

q 2 



l 



K 



(34) 



As a short digression, it is instructive to contrast this result with the analogous expression 
one would obtain for the structure F^ v which leads to an extraction of the magnetic dipole 
moment fi p . The crucial difference is that in this case a nonzero contribution survives in the 
chiral limit. Specifically, a perturbative 1-loop diagram in the background field leads to a 
contribution of the form, 



n 



i 



8tt : 



(e u - e d )F^ In 



-q 2 



+ ■ 



(35) 



Subleading power corrections have been ignored here. However, it turns out that such 
contributions generically have a form similar to those in (|34"D with coefficients m^d 



m q , 

and therefore vanish in the chiral limit m q — > 0. One may also check that the first subleading 
terms of 0{m° q ) actually vanish identically, and therefore the perturbative piece serves as 
the dominant contribution to fi p . An interesting corollary is that, while certain power 
corrections are closely related as above, there would appear to be no simple proportionality 
relation between the electric and magnetic dipole moments. 
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III. MESONIC SPECTRAL FUNCTION AND CONSTRUCTION OF THE SUM 

RULE 



In order to extract a numerical value for the p + EDM from the OPE, we assume as usual 
that Il + satisfies a dispersion relation (ignoring subtractions) of the form, 



OC 



n( 8 ") = -/ (36) 



n Jo (a — q 2 ) 

which we then saturate with physical mesonic states (p + , and excited states with the same 
quantum numbers which we denote collectively as p'). To suppress the contribution of 
excited states, we apply a Borel transform to n + , which we define, following [HJTH], as 



MI+ ee lim s ,_ 00 , s/n=M2? -— — -- n+( S ) = — j o dae-°' M ImIT», (37) 



where s = —q 2 . 

The phenomenological side of the sum rule may be parametrised by considering the form- 
factor Lagrangian which encodes the effective CP violating vertices (see Fig. 2). This has 
the form C = J2 n fnS(q)O n S(q), where f n is the form factor, S(q) is the on-shell propagator 
for p + or one of its excited states, and O n is the operator corresponding to the induced 
vertex. 

As mentioned above, there are a priori two such operators which need to be considered 
at lowest order, F^ u q 2 and q a q[^F u ^ a . As noted above, the second structure vanishes on-shell 
and thus does not enter the form-factor Lagrangian. Another on-shell T-odd form factor, the 
magnetic quadrupole moment, would appear only at the next order in momentum transfer, 
i.e. in front of the structure proportional to d\F^ v |15 |. Since we work only to linear order 



in photon momentum, the magnetic quadrupole moment cannot be recovered from the OPE 
form (|34|), and so we omit it on the phenomenological side as well. 



F 



p' 

FIG. 2. Mesonic contributions to the current correlator in an external electromagnetic field. 
Possible excited states with the p + quantum numbers are denoted generically by p' . 

Consequently, for comparison with the OPE, we have on the phenomenological side in mo- 
mentum space, 

U+jpto.) = 2 f(q 2 )F flu + ■■■, (38) 

where, since we work outside the dispersion relation we may add polynomials in q 2 to ensure 
transversality in the chiral limit, and optimum behaviour for large q 2 , without affecting the 
physical spectral function ImU. We then find that the function f(q 2 ) takes the form, 
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(q 2 -m 2 ) 2 ^ (q 2 - m 2 p )(q 2 - m 2 n ) ^{q 2 - m 2 n )(q 2 - m 2 m )' 

In this expression, A is the dimension 2 coupling defined in terms of the transition amplitude 
for the vector current to go into the p + state, (0|j M |p + ) = XV^, where V p is an appropriate 
vector, while fi is associated with the p + EDM, and f n and f nm correspond respectively to 
transitions between p + and excited states, and between the excited states themselves. 
After performing a Borel transform on /, we obtain^ 

M 4 M z „ m„ — rrip 

+ E -r-2 kr-2 - e -^ +m ^ M2 . (40) 

^ (m 2 n - m 2 p ){m 2 m - m 2 ) 

Since the gap from m p ~ 0.77GeV, to the first excited state m p i ~ 1.7GeV is large, we shall, 
as in [18], ignore the continuum contribution as it is exponentially suppressed. Thus we may 
write, 

= {if + w) e -"' IM < 41 > 

where A is an effective constant of dimension 2. 

We are now in a position to write down the sum rule for the coefficient of F pu . From the 



Borel transform of ([34]), and also (|38"JiT), we have 



A 2 /i + AM 2 = \mj{e u - e d )M 4 e^ 2 <0|gg|0> - ± U + K _ ^ + . . . ( 42 ) 

This is our final result for the CP-odd sum rule, and will be investigated numerically in the 
next section. 



IV. NUMERICAL ANALYSIS 

The coupling A present in (f42"D may be obtained from the well known mass sum rule in 
the CP even sector. In this case, there is no need to consider a background electromagnetic 
field, and the sum rule takes the form (see e.g. 

A ~m 2 /M 2 _ _ 



g pi 



m 2 M 2 4vr 2 







e 


s /M 2 


1 


h 




7T 2 / 


4tt 2 


71 


+ 3M 4 I 



+— ( 1 + - + T^FI ( (0|^G 2 |0) + 24{0\m q qq\0) ) ) . ( 13) 



3 An alternative derivation of the Borel transform in this context, using double dispersion relations 
in order to parametrise I mil, was presented in [11|. 
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Since there is no background field, the leading term is a single pole contribution, and we 
include a continuum term shifted to the right-hand side starting from the p' threshold at 
s ~ 1.7GeV. Note also that A, as defined above, is related to the dimensionless coupling g p 
associated with the width of the resonance, via A = m 2 p /g p , so that g p is dimensionless. 

The physical EDM parameter d p may be obtained by normalising the form factor fx, 
introduced above, by the p + mass. Furthermore, it will be convenient in what follows to 
define an additional parameter d via the relation, 

d A = d^9(e u -e d ). (44) 
m p m p 



We shall now study the sum rule (|4*2|), making use of flUf ) to remove the A-dependence. 
Its helpful to consider the various contributions to (f42|) in turn. At the most naive level, we 
can ignore the 0(1/M 4 ) corrections in (f4"2"D, and also the continuum in (f43|). Taking ratios 
one finds, 

7 2vr 2 x(0\qq\0) 

X ~ m 2 (l + « s /7r + 7r 2 (0 4 )/M 4 )' 1 ' 

where we have defined (O4) = (0|^G 2 |0) + 24(0|m 9 gg|0) for convenience. For numerical 
calculation we make use of the following parameter values: For the quark condensate, we 
have 

(0\qq\0) = -(0.225 GeV) 3 , (46) 
while for the condensate susceptibilities, we have the values calculated in |I7|] and |jI8fl , 



X = -5.7 ±0.6 GeV" 2 [17] (47) 
K = -0.34 ±0.1 @ (48) 



£ = -0.74 ± 0.2 (49) 



Note that Xi which enters at 0(1/M 2 ), since it is dimensionful, is numerically significantly 
larger than k and £. With these parameters, the result for d\ is shown in Fig. 3, where we 
have also used the 1-loop running coupling a s (M) with two flavours, normalised to 0.34 at 
M T . Note that the stability at large M 2 is an artefact of the cancelation of the leading M 
dependence in (f4~5D. One should expect this relation to have reasonable accuracy only in the 
range M 2 ~ m 2 . 

To observe the effect of the 0(1/M 4 ) corrections, we now need to address the issue of the 
unknown constant A in (|4"2"D. Relative to fi, there is a suppression factor of M 2 /(sq — m 2 ) 
associated with A, which near M 2 = m 2 is ~ 0.25. Although this is to be summed over all 
the excited states, we shall use this as justification to treat A perturbatively, and solve for 
it in terms of f\ using (^), but ignoring 0(1/M 4 ) terms. Its convenient to do this via first 
pre-multiplying (f42D by M 2 and then differentiating by 1/M 2 . One obtains the relation 
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FIG. 3. The p + EDM parameter d as a function of M 2 according to various components of the 
sum rules ( |45| ) and ([5l|). 

which vanishes when M 2 = m 2 . Substituting this back into fl42|), we can isolate d by 
taking the ratio with the quantity obtained by pre-multiplying ( f4"3"D by e s °^ M2 , and then 
differentiating by 1/M 2 , 

d = 2t2 (0m0}(s + M 2 -ml) (jL-±( 1 + K -i)+..\ (51) 

(s Q (l + a s /7r) + tt 2 ( So + 2M 2 )((9 4 )/M 4 ) l,M 2 M 4 y 4^ J ' 1 J 

To study the various contributions to (|5"T|), we first set all the 0(1/M 4 ) corrections zero, 
and the result, denoted d 2 , is shown in Fig. 3. We see the expected 1/M 2 behaviour so that 
there is no stability region, although the relation is, as one would expect, very close to di 
(H) near M 2 = m 2 . 

The leading correction at 0(1/M 4 ) may be isolated by setting n = £ = in flolf) . The 
result, ofa, is shown in Fig. 3. The presence of the 1/M 4 term induces a transition region in 
the M 2 dependence. Note that the numerical similarity with d 2 over the relevant range of 
the Borel parameter M 2 is in part due to the compensating effect of the 0(1/M 4 ) correction 
in the denominator of (|51~D . 

Finally, we can obtain an estimate of the corrections associated with k and £, by plotting 



the full expression in (|5T~D , which is displayed in Fig. 3 as d 4 . We see that including these 
corrections has little effect on the behaviour or stability of the sum rule. This is encouraging, 
as the precise numerical values for k and £ are uncertain to a larger degree than that of x- 
Extracting a numerical estimate for d, and an approximate error, from the stability region 
M 2 ~ 0.3 - 0.8GeV 2 in Fig. 3, we find the result 

d " = (2M °- 8H (lG#' (52) 

for the EDM of p + , where e = e u — Cd is the positron charge. As is clear from Fig. 3, the 
dominant contribution arises from the term proportional to the susceptibility \i an d thus 
the result is essentially linearly dependent on the value (and error) for this parameter. 
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It is interesting to note that comparison with the naive estimate (||) implies a value for 
the effective scale A of order 1 GeV. This is very close to a similarly defined scale which 
would effectively appear in the chiral logarithm estimates for d n , Eq. (Q). In this sense, we 
can conclude that our result is in the expected range. One final point to note is that if we 
return to Eq. (|5l|), which is written in terms of the condensates, and re-express the final 
answer in slightly different units, we can directly relate the EDM to the vacuum topological 
susceptibility correlator K as calculated in 

d p = 2.2 x lO- 3 e , ffii,., rn u m d = % 1Q _ 3 Q (53) 

p (lOOMeV) 5 (m u + m d ) 2 (lOOMeV) 5 v ' 

Note that we have used a normalisation (lOOMeV) which is adapted to the small size of 
m*, and it is this which accounts for the small overall coefficient. This factor is essentially 
hidden in the result presented in ([52]). 



V. DISCUSSION 

Throughout the calculation we have intentionally kept m u = m d , knowing that the 
correct mass dependence is m*. It is easy to see, however, that if m u ^ m d the calculation 
does not automatically restore m* since, for example, the up-quark bilinear combination in 
Eq. (p~3|) comes with a coefficient proportional to m d . In the final result it would induce a 
contribution which would not vanish in the limit m u — > 0. This means that if m u ^ m d there 
should be additional contributions which would combine with the rest to form an overall 
m^-dependent result. 

At the same time, we recall that one can use a chiral transformation in the QCD La- 
grangian and rotate the ^-parameter to stand in front of the quark singlet combination 
m*(m'75W + di^d). It is clear that in this situation the ^-dependence for any physical ob- 
servable will arise together with the correct mass dependence and will disappear at m u = 0. 
The answer to this "puzzle" lies in the chirally non-invariant form of the quark current which 
we associate with p + . Using this form of the current, additional contributions must arise 
which are associated with vector-axial-vector current mixing^. In other words, the purely 
vectorial current is not diagonal due to the chiral anomaly, and one needs to consider all the 
(jvjv)i UvJa), and (jaJa) correlators to obtain a well-defined projection onto p. However, 
a more elegant approach would appear to involve a direct diagonalisation of the current. To 
see how this might be achieved, let us write down two forms of the QCD Lagrangian with 
an external vector source coupled to isovector quark current: 

d = m u uu - m d dd + 9-^G a G + V^u^d + V*d^^u (54) 

£ 2 — . . . _ m u uu — m d dd — Om^ui^u + dijsd) + V^u^^d + V*d^^u (55) 

where the ellipses stand for the standard kinetic terms for the gauge and quark fields. In 
the absence of the external current C± and £2 are equivalent (we consider 9 to be small 



4 We thank A. Vainshtein for discussions on this point. 
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and work only to linear order) . The presence of the external current in the form written in 
the Eqs. ([54]) and (|55"D makes £\ and £2 explicitly inequivalent. The same chiral rotation 
transforms £\ to £2 plus an extra term 

£1 — > £ 2 + iQ d , u iy iJ ,vn^d - V*d^ 5 u). (56) 
m u + m d 

In the limit of m u = 0, this extra term contains 9 explicitly which will then enter in the 
physical amplitudes bilinear in V. Thus we need to bear in mind that in the presence of 9 
the choice of the current for further use in QCD sum rules is not "automatic" , if one wishes 
to avoid mixing with other contributions. In the case of p + with m u 7^ m^, the G?7 M w-current 
should be used only in the basis where 9 is completely rotated to the quark mass term. In 
the basis where 9 enters in front of GG, the current includes additional axial-vector pieces 
which restore the correct quark mass dependence in the final answer. As an additional check, 
one can calculate the next order ~ # 2 -corrections to CP-even observables and observe the 
dependence of the result on the choice of the current. 

Of course, the calculation of the electric dipole moment of p does not have direct ex- 
perimental implication. Our main motivation for calculating this quantity was to test the 
possibility of applying the QCD sum rule approach to the problem of EDM(#). We would 
like to mention here that the idea to consider the EDM of the p meson resulting from the 



EDM of quarks was used in [19 



Returning to the problem of d n (9), it seems clear that there are a number of additional 
difficulties which may be encountered. One of them refers to the correct choice of the nucleon 
current at 9 7^ 0, as this choice can be ambiguous even in the normal CP-conserving case 
|20| . Another difficulty is related to the necessity for a simultaneous treatment of the mass 



operator and the electromagnetic form factors. This is because in the presence of 9 the mass 
operator develops an imaginary part which can influence the answer for d n (9). Nevertheless, 
this calculation appears feasible and work in this direction is currently in progress. Only 
then can one have a reliable means to interpret d n directly in terms of the high-energy 
parameters (CP-violating phases in the soft-breaking sector and masses of superpartners in 
the case of the MSSM). 

In conclusion, we have demonstrated that QCD sum rules can be used for the calculation 
of CP-odd electromagnetic form factors induced by the theta term. The result for the EDM 
of the vector meson, calculated in this way, is stable and numerically dominated by the 
vacuum magnetic susceptibility. The set of correlators which appear in the OPE part of the 
QCD sum rule were calculated via the use of the anomaly equation, in a similar manner 
to the calculation of topological susceptibility. In this way, a direct relation to the U(l) 
problem becomes apparent as the total result vanishes in the limit m n — > m n . 
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